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Quantum linear response theory considers only the response of a closed quantum system to a 
perturbation up to first order in the perturbation. This theory breaks down when the system subjects 
to environments and the response up to second order in perturbation is not negligible. In this paper, 
we develop a quantum nonlinear response theory for open systems. We first formulate this theory in 
terms of general susceptibility, then apply it to deriving the Hall conductance for the open system 
at finite temperature. Taking the two-band model as an example, we derive the Hall conductance 
for the two-band model. We calculate the Hall conductance for a two-dimensional ferromagnetic 
electron gas and a two-dimensional lattice model via different expressions for d a (p), a = x,y,z. 

The results show that the transition points of topological phase almost remain unchanged in the 
presence of environments. 
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I. INTRODUCTION 

Applying a small perturbation to a quantum system, 
we may compute the response of some operators of the 
system to the perturbation, for example the electric cur¬ 
rent density as a response to an electric field. When 
consider the response only to first order in the perturba¬ 
tion, the theory is the so called quantum linear response 
theory. 

The Kubo Formula is an equation which expresses the 
linear response of an observable quantity due to a time- 
dependent perturbation, which has been widely used in 
condensed matter physics EH3 since it was first derived 
by Kubo in 1957 [|[. However, Kubo’s theory is only valid 
for an equilibrium closed system in the linear regime @. 
In recent years, linear response theory for open systems 
has attracted more and more attentions in biophysics, 
nano-physics and condensed matter physics. A linear re¬ 
sponse theory based on the master equation and 

the hierarchical equation of motion 13|, [14] has been de¬ 
veloped for open systems. The key issue and difference 
of those approaches are how to get the reduced density 
matrix of the open system—while the first approach ob¬ 
tains the density matrix by master equations, the second 
by the Hierarchical equation. Both approaches need to 
calculate the density matrix or the dynamics of the den¬ 
sity matrix. This means we have to trace out the envi¬ 
ronment first, then calculate the response—a treatment 
on nonequal footing for the system and environment. It 
is worth addressing that, based on the reduced density 
matrix, a response theory for the open system was devel¬ 
oped, we refer the readers to EE] for details. 

Treating the system and the environment on equal 
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footing, we here develop a quantum response theory for 
open system by a different approach. This approach does 
not need to get the reduced density matrix and the the¬ 
ory is beyond the linear response regime EMI. Our 
treatment is not limited to a specific system. Rather, 
it is applicable to any strength of perturbation and open 
quantum systems. Hence it is suitable for strong external 
fields applied to topological insulators El- We will ap¬ 
ply this theory to topological insulators (TIs) described 
by the two-band model and study the effect of environ¬ 
ment on the Hall conductance at both zero and finite 
temperatures. 

Topological insulators are a broad class of unconven¬ 
tional materials that are insulating in the interior but 
conduct along the edges [Tt 1 - [20I . Over the last decades, 
topological insulators have attracted a great deal of inter¬ 
est due to their interesting features and possible applica¬ 
tion in quantum computation. Recently, efforts have also 
been made to investigate topological insulator for open 
systems, e.jr density-matrix Chern insulators by ther¬ 
mal noise |2ll . [2^ |. zero-temperature Hall conductance 
subjected to decoherence EEL topological order by dis¬ 
sipation [23}. These stimulate us to develop a response 
theory for open systems to high order in perturbation, 
and apply it to study topological insulators subjected to 
environments. 

The remainder of the paper is organized as follows. In 
Sec. II, we extend the linear response theory of closed 
system to open systems and nonlinear case. In Sec. Ill, 
as an application of our theory, we derive the finite- 
temperature Hall conductance for open systems and ex¬ 
emplify it into the two-band model. The Hall conduc¬ 
tance for a two-dimensional ferromagnetic electron gas 
and a two-dimensional lattice model is also discussed in 
this section. Discussion and conclusions are given in Sec. 
IV. 
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II. NONLINEAR QUANTUM RESPONSE 
THEORY FOR QUANTUM OPEN SYSTEMS 


Consider a quantum system described by a time- 
dependent Hamiltonian Hs(t) coupling to environment 
Hr in an external field H e (t), the total Hamiltonian 
Hx{t) reads, 

H T (t)=Hs(t) + H R + H S R+eH e (t), (1) 


where H$r denotes the interaction between the system 
and the environment, e stands for the field strength, and 
H e {t ) = — Y! fv{t)C v - To simplify the notations in later 

V 

discussions, we set H{t) = Ho(t) + HsR , Ho(t) = Hs(t) + 
Hr, and use the notation, Lx(t)p(t) = [Hx(t),p(t)] with 
X = T,e, 0, or left blank. We are interested in the re¬ 
sponse of the open system to the external fields /„, which 
will be derived in the following. 

The total density matrix px(t) satisfies the Liou- 
ville equation pr(t) = — ^-Lt(£)pt(V) = — jrL(t)pr(t) — 
ieL e (t)pT(t). By dividing the total density matrix into 
two parts p(t) and pTe(t) with pit) satisfying p(t) = 
— j^L(t)p(t), we have 


Pe(t) = --^eTtr 


g(t, u)L e (u)p T (u)du, 


( 2 ) 


where p e (t) was defined as p e (t) = Tr R {pTe{t)) and 
PTeit) is defined as a change of the density matrix due 
to the external field. Trx denotes a trace over X. 

Eq. j2j) completely describes the influence of the ex¬ 
ternal fields on the system subjected to an environ¬ 
ment. g(t, u) is the chronological time-ordering operator, 


g{t,u) = T + exp -j-f^dsL(s) , (t > u). Now we can 


define the general nonlinear response tensor Xnv{^i u ) f° r 
the open system via the change of the expectation value 
of a system operator F ^ caused by the external field, 
(F„(t)) e = Tr s (Pe(t)Fj, 


( F »(t)) e = '52 [ duXnv(t,u)f„(u), 

,, Jtn 


(3) 


whole system is in an equilibrium state p eq at tempera¬ 
ture T. The main task of the general nonlinear response 
theory is to calculate the susceptibility, 

pOO 

X/^M = / dte lult x^( 0 ,-t) 

Jo 

poo 

= l -ej dte iut TX S R{F^e-^ Lt [C v ,x{-t)}} 

( 6 ) 

with an assumption that H e is turned on at i 0 —> —oo and 
the steady response outputs at t = 0. Defining q v {t) = 
Tr Re~^ Lt [C u ,x(—t)\, we rewrite the susceptibility as, 

i r°° 

X/^M = -^e J dte lLUt Ti s [F^q^ft)]. ( 7 ) 

By the modified Laplace transformation g(ui) = 
/ 0 °° dte lut g(t) (HI, we find, 

OO 

X^M = eyW(w) (8) 

n —2 

where x${u) = j-Ti s [F^qi n) (w)] with 

pOO 

<li n) (u) = / dtexp(iut)ql n \t), 

Jo 

q^Ht) = TV R e~^ Lt [C v ,p eg ], (9) 

and 

ql n) (t) = (-l)" +1 Tr R {e-i Lt [Cy, [ ‘ dti ■ • • dt n 

J —oo J —oo 

•s( — t,tl) ■ ■ ■ S(t n -2,t n -l)peq]}, {n > 2). (10) 


Here the first term in Eq. (|5]) represents the linear re¬ 
sponse, while the others are nonlinear responses of the 
open system. Eq. (0 suggests that in order to get 
the nonlinear susceptibility, we have to calculate q v (ui) 
through q v {t). The time evolution of q v {t) is given by 
(see Appendix B), 


where the nonlinear response tensor X/iv(t, u) is given by 

u ) = t^T rsR{F^,g{t, u)[C v , z(u)]}, (4) 

with 

00 pU ptn — 1 

x[u) = p(u) + ^2 ( -£ ) n / d<l / dtn 
n=1 ■'*0 -' t o 

■s(u,ti) ■ ■ ■ s(t n -i,t n )p(t n ), (5) 

and s(t,u) = j-g(t,u)L e (u). 

In the following, we will restrict ourself to consider 
a time-independent system. The formalism can be eas¬ 
ily generalized to time-dependent systems. Suppose the 


i 

q„{t) =--[H s ,q v {t)\+ I c(t-T)q v (r)dT + K{t), 

(U) 

where the kernel c(t ) and term K(t) are given by c(t) = 
—Tr R[Le~TiQ Lt QLpR/h 2 } and 

K{t) = -^Tr R [Le- iQLt / h QA(0)\+Tr R e-* Lt A{t) 

——Ttr[L f dTe~^ < ^ L< ' t ~ T ^Qe~T [LT A(r)], 

™ Jo 

( 12 ) 

respectively. A (t) = [C v ,x{—t)\, and x(t) is defined by 

Eq. 10. 
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In Appendix C, we use an example—a single-mode cav¬ 
ity system coupled to environment in a time-dependent 
external field—to illustrate the nonlinear response when 
the external field is not weak. In addition, we show that it 
is difficult to derive analytically the nonlinear response of 
a general open system to a perturbation at finite temper¬ 
atures. Considering the fact that two-dimensional topo¬ 
logical insulators described by the two-band model has 
been experimentally observed jl7l420|. we will focus on 
the linear response theory for an open system in the fol¬ 
lowing. 

We assume the system and environment be ini¬ 
tially in their thermal equilibrium, p( 0) = p eq = 

e -^H/Tr SR e~P H . To calculate the linear response, we 
take the first term in Eq. ([5]) as x(u) = p eq ■ Further¬ 
more, we take p eq = ps ® Pr + O(Hrr) for simplicity, 
the term K(t ) in Eq. (fill) vanishes since Q[C„. p eq ) = 
Q[C Vl Ps)pr = 0. After the modified Laplace transfor¬ 
mation for Eq. CD and expanding it to second order in 
the system-environment couplings Hsr , we have 

- -^Lsqv(uj) + [iw + c(u)]q v (ui) = -q v { 0), (13) 


where 


c w = 


-1 


dt{L SR H I SR {t)) R exp[i(w - L s /h)t], 


(14) 


with T-Lg^t) = exp[—iLot/h\HsR . Here and hereafter, 
the perturbation parameter e is absorbed in C v . Eq. m 
and Eq. (0 is one of the main results in this paper. Eq. 
0 is the response of the open system to the external 
field, and Eq. (fT3l) is the key element to calculate the 
linear response. 


III. APPLICATION TO HALL CONDUCTANCE 
IN THE TWO-BAND MODEL. 

We now apply the quantum response theory for open 
systems to a two-band model. The model Hamiltonian 
is, 


Hs = 'Y^H s (p),H s (p) = e(p) + ^2 d a (p)a a , (15) 

p a.=x,y : z 

where e(p) = p 2 /2m* denotes the kinetic energy with 
the band electron effective mass m* and u a are the Pauli 
matrices, p = (p x , p y ) stands for the Bloch wave vector 
of the electron. 

The two-band system is an idealization but can be re¬ 
alized approximatively with ultra-cold atoms 25 using 
different techniques as, e.g., super-lattices [26 b 2£ with 
hole-related band suggested by Raghu et al. in 29] and 
in spin-1/2 electrons with the spin-orbit coupling [3j. Be¬ 
sides the possibility of experimental realization, the two- 
band model is also interesting as a simple model in con¬ 
densed matter physics. 


The eigenenergies of the two-band Hamiltonian are, 
E m (p) = e(p) + md(p), with m = ± and d = 



where cos 9 = d z /d, tan (j> = d y /d x were defined. 

The Hall conductivity tensor a^ can be calculated 
through the current density J ;/ , (t) in the /x-direction 
(/x = x,y,z ), as a response to the external electric field 
E v (t) = RefEy exp(xwt)] in that direction, the current 
density reads, 

J y (t) = lim Re{a^ v (uj)E v exp(xwt)]. (17) 

oj —>-0 


Consider a total Hamiltonian in the momentum space 


H T (t) = Hs + H R + H S R + H e (t), (18) 


with Hr = 'Y/,hbjjb^bj 1 HsR = + H.c., and 

j j 

H e (t) = — Y)PvE v (t). Here Hs is given by Eq. (fl5)l . 

Eq. (1181) denotes the Hamiltonian of a system composed 
of many particles with charges q :j and position operators 
r 5 = 'ihgyj subjected to environmental noises Hr, Py is 

the polarization Py = ifylj jtt ■ I n the linear response 
j Pu 

theory, the Hall conductance takes [30j 


ff M y(o;) 


1 r°° 

=- h J Q dt([P u (0),J,(t)])e^ 
J^° dt([M0),Mt)])e iut , 


(19) 


where J u = P u /V and V denotes the mode volume for 
system. Using the relations J = ev/V and n M (f) = 
e lHt v ll e~ lHt , we obtain, 




-e 2 r°° 

=j~vj o 

p 2 l-oo 

=T^V l (20) 

+ O(Hsr) 


huV 


■Tr s [v^q v (u)]. 


Here we set C v = v v , therefore q v ( 0) = [v v , Ps\ + 0(Hsr) 
. Noticing Eq. (fl8l) and taking the Lorentzian spectrum 
density J(u>) = ^ (31] into consideration, under 

the Markovian approximation we obtain from Eq. m 


c(u)q v (u) = r[2cr_gy(w)cr + - cr + <T_q„(uj) - q v (u)<r + a_]. 

( 21 ) 
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Substituting Eq. m into Eq. m and inserting the 
identity \m(p)) (m(p) \ = I into the result, we obtain 

m 


0 — Sis,nm H - ^Qi/,nm\P^ ^nip) ■^m{.P)\ 

+ 2TH ^ ' { a —) n iqv,ij(< J +)jm ~ Ff,nm ~~ F'v^mni 

P,ij 

( 22 ) 

where the coefficient F v ^ nm = hT^2q Vtn j(a+a-)j m , 

p,j 

Sv,nm = fi[fm{p) ~ fn{p)\v„ t nm with the definition y mn = 
(m(p)\y \n(pj). Here we have applied the fact that ps = 

E 1 3ip)) fj(p){j(p)\, and fj(p) = l/{exp \P(Ej(p) - p)] + 
j 

1} is the Fermi-Dirac distribution function with p the 
chemical potential, /3 = 1 /ksT, ks the Boltzmann con¬ 
stant, and T the temperature. Tedious but straightfor¬ 
ward algebra yields (m yf n), 


$is,mnA nm HC B nm COs26 Si/, nrn (hu -f- iAr)IzT 

'-fi/ nm — o q 

2[D nrn 2(huj) + 8i(Hcu) hT — ihT 6^00820] 

(23) 

where the coefficients take e nm = E n ( p) — E m ( p ), A nm = 
4i(fiw) 2 - tuo (llftr + 4ie„ m ) - M\7ifir - 6e„ m ), B nm = 
k>u.rirn.(Ed T iST) T S u , mn {huj -(- ifiT 2e nm ), and D nrn — 
-2hu[S(hT) 2 + e 2 J - ihmm 2 + 3e 2 „J. 

The diagonal elements of q u {^) is not listed here, since 
it has no contribution to the Hall conductivity. In the 
weak dissipation limit, T —> 0, we can expand q u ^ n m in 
powers of T. To first order in T, 5 „ inm can be written as, 

Qv,nm = Qv,nm 4“ (24) 

where the zeroth- and first-order of <?„ jn7n take, 


9 

9 


(o) 

iy,nm 


(i) 

v,nm 


iS v 


klO Cnm 

A(hio - e nm f 


h{6)S„, r 


4(SV-eL)’ 


(25) 


respectively. g{6) = 5 — cos 2 0,h(9) = 1 — cos 26. Sub¬ 
stituting Eq. (HH1) into the third equation of Eq. CU1) , we 
have 


M = cr^J (w) + nr crW (w), (26) 


with the zeroth- and first-order Hall conductivity at 
finite-temperature T, 



a 


(i) 

fllS 


M 



HujV 

p,m^pn 


(m{p)\vy, | n{p}) 


y{0)S un m 
A{tno - e nm ) 2 


+ 


h^O^Svmn 
A(h 2 u 2 - e 2 nm ) 


(27) 

Finally, we can obtain the finite-temperature Hall con¬ 
ductance for the open system in the weak dissipation 


limit, (i.e., T —> 0, see Appendix D) 


+ cr$ = C-, 


(28) 


where 


(0) e2 [ /A 
O'E = 2fr T ( fc ) £ «/37 77^2 ; 


(2?r) 


^ I T ® (l 9(0) 

ol,P ,‘7 


a/3 


Cl CL -y d'y 


d 2 k 


(2vr) 2 
(29) 

with t(LA — l/+( fc )-/-( fc ) 1 dd a ddp , n 

Wlttl T{h) — dk„ dk v a 7; ~ 

lm[(+(p)\a a \-{p}) (+{p)\ap |-(p))]. C in Eq. (EED de¬ 
fines the Chern number of the open system. Here we 
have used )T) —>• ^ h yi f dp x dp y and p = hk. Note that 

p 

the correction in general is complex, it contains a real 
part (the first term) and an imaginary part (the second 
term). In the following, we will present two examples. In 
the first, the correction to the Hall conductance from the 
system-environment coupling is real, while it is imaginary 
in the second. The two examples together exemplify the 
effect of the environment on the Hall conductance. 
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FIG. 1: (Color online) The zero-order conductivity ai°J [(a) 
and (b)]( i.e., F = 0) and the total conductivity a xy [(c) 
and (d)] as a function of /3o (meV-nm/ft). The zero-order 
conductivity is exactly the conductivity of the closed system. 
It is worth addressing that the imaginary part of the Hall 
conductivity is zero in this case. Parameters chosen are T=0.2 
meV/ft, Ao = 5 meV-nm/fi, p = 1 meV, ho = 2 meV, in * = 
0.9m e , m e is the mass of electron. The temperature T = OK 
for (a) and (c), and T = 300A' for (b) and (d). 


Example 1 .— Consider a two-dimensional ferromag¬ 
netic electron gas in the presence of both Rashba and 
Dresselhaus spin-orbit couplings m, m- This sys¬ 
tem can be described by Hamiltonian Eq. m with 
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d x = Xohky — (3ohk x , d y = —Ao hk x — f3ohk y , and d z = ho. 

By the use of Eq. (T251> , we numerically calculate the Hall 
conductance and plot the results as a function of /3o in 
Fig. ED Fig. |T] (a) shows a phase transition at /3o = ±Ao. 
When /3 q > Ag, the Chern number of the closed system 
at zero-temperature is 0.5, while for /3g < Ag, the Chern 
number is —0.5. This is in agreement with the analytical 
results of Hall conductivity of the closed system given by 
axy 1 (T = 0) = isgn(/3g — Ag) in units of For the open 
system, the plot shows that the phase transition can still 
emerge. This can be found in Fig. [T] (c), the Chern num¬ 
ber at zero-temperature is about 0.5 when fig > Ag, while 
for 0g < Ag, the Chern number is about —0.5. A critical 
value /3o = ±Ao at which the Hall conductance changes 
abruptly can be found. This suggests that the topolog¬ 
ical phase transition survives for open system. We also 
find that at the critical point Ag, a xy increases compared 

to <7®, while it decreases at — Ag. Except the two critical 
points. The effect of the environment on the conductance 
is almost zero. Thermal fluctuation diminishes the dif¬ 
ference of the Hall conductivity between different phases, 
but it does not change the nature of the phase, see Fig. |T| 
(b) and (d). It is worth addressing that the imaginary 
part of a2y vanishes for in this model, this is different 
from the example below. 



FIG. 3: (Color online) The Hall conductivity of the closed 
system [(a) and (b)] and the imaginary part of the first 
order Hall conductivity Imfcr^'D] [(c) and (d)] of the open sys¬ 
tem. We plot them as a function of t a (meV) and <5 (meV) 
at finite temperatures. We address that the real part Re[cr)u)] 
of the Hall conductivity of the open system is zero in this 
case. Parameters chosen are F = 0.1 mev/fi, u> = 0.2 mev/fi, 
p = 1,5 = 4, l = 1, mo = mi, T = 30A' for (a) and (c). 
T = 300A for (b) and (d). 



FIG. 2: (Color online) The zero-order Hall conductivity ai°J 
[(a) and (b)] and the imaginary part Im [cr$] [(c) and (d)] of 
the first-order Hall conductivity as a function of t a (meV) and 
5 (meV) at zero temperature. Note that the real part Re[cr£]?] 
of the Hall conductivity of the open system is zero in this case. 
Define two coefficients, mi = 4m-(-1 or 4ru-|-4 (m=0,l,2, •••), 
and m 2 = 4m + 2 or 4m + 3 (m=0,1,2, • • ■). Parameters 
chosen are T = 0.1 mev/fi, a; = 0.2 mev/fi, p = 1, q = 4, l = 1. 
mo = mi, for (a) and (c). mo = m 2 , for (b) and (d). 

Example 2. — As the second example, we consider 
tight-binding electrons in a two-dimensional lattice de¬ 
scribed by the Hamiltonian 


H = 


(*.j) 


„f„ e iSi 


„ctcje*®* 

<*»j> 


where Cj is the fermion annihilation operator on the lat¬ 
tice site j, t a and tb denote the hopping amplitudes along 
the x — and y — direction, respectively. Consider two 
branches coupled by |Z|—th order perturbation, the effec¬ 
tive Hamiltonian then takes the form Eq. (1151) with d x = 
6cos(k y l ), d y = <5sin (k y l), and d z = 2t a cos(k x + 2n^mo), 
where p and q are integers. 5 is proportional to (is the 
order of) fj^. 

In Fig. OH we show numerically the Hall conductance 
at zero-temperature for the open system as a function of 
t a and <5, the Fermi energy is set in the gap. We find 
that the Hall conductance a^y and Im [a^J ] change its 
sign when t a crosses zero. The topological phase tran¬ 
sition at zero temperature survives in the open system, 
this can be observed by examining the Hall conductance, 
which changes from —0.5 and —0.1 to 0.5 and 0.1 when 
mo = mi [in units of e 2 /fi, see Figs. [2] (a) and (c)]. 
Similar changes from 0.5 and 0.1 to —0.5 and —0.1 for 
mo = m 2 [see Figs. [2] (b) and (d)] are found. We here 
address that the behavior of the imaginary part Im[crp'j] 
of the Hall conductivity for the open system is the same 
as afnl for the closed system [see Figs. 03] except their 
amplitudes. This means that the topological properties 
of open and closed system are same. This observation 
can be explained as follows. We analytically calculate 
the Hall conductance at zero temperature and obtain, 






(30) 


(31) 
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where A = -Jfe,B = ni3 Z + M 2) ’ M = < + S 2 >0. 
Especially, we find 


T (°) = 
' fiiy 

tW = 


^ 2^ Sgn(£ a )sgR[$ momi ^mom2]) 


i3zr 

I2u7 c 


( 0 ) 


(32) 


when <5 = 0. This gives the phase transition points for 
the closed and open system at zero temperature with 
(5 = 0 [see Fig. [2] (a) and (c), (b) and (d)]. In contrast 
to the system in the last example, the correction to the 
conductance due to the system-environment coupling is 
imaginary, which wc will refer to environment-induced re¬ 
actance. This environment induced reactance describes 
the energy exchange between the system and the environ¬ 
ment, reminiscent of the non-Markovian effect. More¬ 
over, from Eq. m we find the reactance Im [cr$] and 
a^°J have the same sign because B > 0. At finite temper¬ 
atures, for example, T = 30K[(a) and (c)] and 300K[(b) 
and (d)](see Fig. [3]), the thermal effect diminishes the 
amplitude of the Hall conductance compared with the 
value at T = OK (Fig. [2]), but it does not change the 
topological phase. 


IV. DISCUSSION 

We have developed a quantum response theory for 
open systems beyond the linear regime. A general non¬ 
linear susceptibility for the open system is derived. This 
theory provides us with a formalism to extend the notion 
of the finite-temperature Hall conductance form closed to 
open systems. This comes into play when studying deco¬ 
herence effects on the Hall conductivity in open systems. 
We exemplify the theory in a two-band model that de¬ 
scribes topological insulators, the results show that the 
environment affect slightly the topological phase tran¬ 
sition at finite-temperature. Although the analysis has 
been restricted to two-band models, we expect that the 
general response theory could be extended to higher di¬ 
mensions and nonlinear Hall conductance for many sorts 
of topological insulators. 

The prediction can be observed in the six-terminal de¬ 
vice (for details, see Ref.[36j). The environment can be 
simulated by the use of Btittiker’s virtual probes [[37, 38|. 
As we show, the phase transition remains in the open 
system. Thus, these quantum plateaus are observable in 
the mesoscopic sample even with environmental noise. 
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Appendix A: The validity of the approximation 

q„(Q) ~ \v u ,ps]. 
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FIG. 4: (Color online) The Hall conductance for closed system 
[(a) and (b)] and open system [(c) and (d)[ given by the third 
equation in (120)1. The purpose of this figure is to show the 
difference in caused by different order of q y (0) in Eq. (IA2I) 
by numerically solving Eq. m • The red-solid and blue- 
dashed lines correspond to g„(0) = [v„,ps] + [v„,Pg^ and 
q v { 0) = [u„,ps], respectively. Parameters chosen are u> = 0.2 
mev/fi, A = 25 mev/fi, p = 1, q = 4, L — 1, mo = 1, T = 30 K, 
5 = 0.01, r = 0 for (a), 5 = 2 mev, T = 0 and (b), 5 = 0.01 
mev, r = 0.1 mev/fi for (c), and 5 = 2 mev, T = 0.1 mev/fi 
for (d). 


Here, we examine the validity of the approximating 
<?!/(0) = Tr R [v v ,p eq ] ~ [v v ,p s \ + 0(H S r) in Eq. (0 . 
To this end, we employ the formula with two Hermitian 
operators A and B, 

e A + B = e A e B e~^ A ’ B ^ ■ e ( hi S h -° rder terms for B) 


Under weak system-environment interactions HsRi we 
can expand p eq in powers of Hsr- To second-order 
in Hsr and setting A = —f3H 0 and B = —/ 3Hsr in 
Eq. ED, we obtain 


gp(0) = [ v„,Tr R p eq ] ~ [v v ,p s ] + [w,Ps 2) ], (A2) 


with 


Ps = Ps[ai{Hs<r+H s a-+ cr + H s a-Hs 
-Hscr+a-Hs - a + H s H s <J-) 
+ci2<7+Hscr-+ a 3 cr + <7-], (A3) 


where the coefficients ai = and a 3 = 

TX P . Here we have used the Lorentzian spectral density 


311- 


In order to check the validity of <fi/(0) ~ \v v ,ps\ in 
Eq. (1T31I . we numerically plot the Hall conductance for 
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closed and open system in the tight-binding model (l30l) 
with different <5. In Fig. |ffj we show a comparison of 
the results with two different orders in < 7 ^( 0 ) given in 
Eq. (1A2I) : the simulation is performed for the Hall con¬ 
ductance given by the third equation of Eq. m with 
Eq. (THU) . Blue-dashed line [closed system (a) and (b) and 
open system (c) and (d)] in Fig. @] is for q l ,{ 0) = [v v ,ps], 
which are in good agreement with the results obtained 

/o\ 

with q„(0) = [v V) ps\ + \p v , Ps q- In addition, we find the 
higher-order terms in g„(0), which have no effects on the 
phase transition point and can be ignored with respect 
to the zero-order term in the two-band model. 


Appendix B: The derivation of the equation for g M (t). 


We apply the projection operator method to de¬ 
rive the time derivative of Quit). Note that q v (t) = 
Trne~ii Lt [C v ,x(—t)] is not a density matrix since the 
operator C v does not commute with x(—t) given by 
Eq. flSJ| . To derive the equation, we define superoper¬ 
ators P(...) = pr <g) Trjt(...) and Q = 1 — P, which are 
projection operators satisfying P 2 = P and Q 2 = Q. By 
the standard procedure for deriving a master equation 
in Ref. JTo, !&| 5 q| , with two time-evolution operators de¬ 
fined by (the initial time is to = 0), 


W(t) = Pe-* Lt A(t), 
Y[t) = Qe~^ Lt A(t), 


(Bl) 


Appendix C: Quantum response beyond the linear 
regime. 

As mentioned in the before, the linear response theory 
is not valid when the external field is not weak. Here 
we present an example to show the difference between 
the linear response theory and the nonlinear response 
theory, it also is an illustration for the essential role of the 
nonlinear response theory. We exemplify the difference 
though n e (t), which is defined as the difference in photon 
number with and without external fields inside a cavity. 
The results show that when the coupling between the 
field and system is weak, the linear response is a good 
approximation, otherwise nonlinear response should be 
taken into account. 

Consider a single-mode cavity system with bare fre¬ 
quency wo coupled to a non-Markovian reservoir, driven 
by an external laser with frequency w^. We assume the 
reservoir modeled by a set of harmonic oscillators is at 
a finite (ambient) temperature. In a rotating frame, the 
Hamiltonian of the total system reads, 

H T = H + H e (Cl) 


with 

H =HAa 1 a + ^ hVL k b\bk + ^ (hglab\. + hgkb k a ^, 
k k 

(C2) 

and 


we have 


H e = Ml(a + aJ), (C3) 


W(t) = —jPLW{t) - jPLY(t) + Pe~i Lt A(t), (B2) 


and 


Y(t) = —~QLY(t) - -QLW(t) + Qe~i Lt A{t). (B3) 
n n 

Solving Eq. IB31 we have 

Y(t) =e~i QLt QA(0) - [ e~i QL{t - T) QLW{T)dT 

h Jo 

+ [ e-^ Qi(t - r) Qe-^ iT A(r)dr. 

Jo 

(B4) 

Substituting Eq. (IB41) into Eq. (IB2I) . we obtain 


W(t) = - —PLW(t) - -PLe~i QLt QA(0) - —PL 


o~i QHt-r 


^QLW (t)cIt - l -PL 


e -j-QL(t-T) nfl -i LT 


Qe n T A(r)dr + Pe h 4 A (t). 
J o 

(B5) 

Finally applying P(...) = pn®Tr R (...) into Eq. IB51 we 
arrive at Eq. CD- 


where A = wo — wl, and flk = Wfc — wl. Q is the strength 
of the external field, a is the cavity annihilation operator, 
and b k and g k are the reservoir annihilation operator and 
coupling constant. In the following we will solve the exact 
non-Markovian dynamics in Heisenberg picture. 

Suppose the system and the environment is initially 
uncorrelated—the reservoir modeled by Hamiltonian 
Hr = £ fc fiw*4 6 fc is in a thermal equilibrium state, 
while the system is in a coherent state. The initial state 
of the total system is, 


p h* 11 R 

p T (0) = |a) {a\ <8> p R (0), p R (0) = iYe -fM R > ( C4 ) 

where |a) is a coherent state and /3 = 1/krT. 

Our task is to obtain exactly the response of the system 
to the external field and compare it with the nonlinear 
response. We use the word exactly to denote that the re¬ 
sponse is not a perturbative result- it is exact, including 
all orders in the external field. With the formal solution 
of bk(t) = e^ HTt bk{0)e~^ Hrt , we rewrite the equation 
for a(t), 

d 

—a (t) = - iuj 0 a ( t) ~ J o a (r) / (t - r) dr 
— ib(t) — iLl ( t ), 



where b(t) = Y^9kbk (0) e iulkt . The memory kernel 
k 

f(r) = '£, \9kfe~ iuikT = f dtoJ ( W ) e~^ 

k J 

characterizes the non-Markovian dynamics of the reser¬ 
voir. 

Because of the linearity of Eq. (IC5I) , a (t) can be ex¬ 
pressed as a(t ) = ui(t)a(0) +v±(t), where a (0) and 
bk (0) are the operators at the initial time. Here time- 
dependent coefficient u\ ( t ) and v\ ( t ) can be calculated 
by Eq. (1C5|) , 

( t ) = -iuj 0 ui (t) - J ui (r) f (t — t) dr, (C6) 


d /'* 

—V! (t) = - iuj 0 v i ( t ) - / Vi (r) f (t-r) dr 
dt J 0 

— ib(t) — ifl (t ), 


(C7) 


with initial conditions u\ (0) = 1 and V\ (0) = 0. 
Vi (t) can be given analytically by solving the inho¬ 
mogeneous equation of Eq. (1C71) . it leads to v± (t) = 

- i fo i b ( T ) + n {t)]ui {t - T)dT. 

Taking the state in Eq. (IC4I) into account, we can 
calculate the change of the expectation value of the cav¬ 
ity photon n(t) = Trs[aJ {t)a{t)p( 0)] in relevance to the 
external field (1C3I) as, 


n e (t) = nr(t) — n(t) = 2R e[u*(t)y(t)] + D(t), (C8) 

where y(t) = —ia*Di(t),D(t) = \Di(t)\ 2 , D\{t) = 
Mfo wi(r)dr. 

We assume that the system coupled to a reservoir 

p \2 

has a Lorentzian spectral density J(w) = ^ p . x2 
till l42j | . In order to examine the validity of linear re¬ 
sponse, we plot the response of the average photon num¬ 
ber to external field in three regimes divided by linear, 
two-order nonlinear and exact response in Figs. [5] and 
6. In non-Markovian regime, e.g., A = 0.2T, we can see 
that the results given by Eq. © under the first-order 
approximation (containing only the first term) in Eq © 
are in good agreement with those obtained by the exact 
response (IC8I) when the interaction strength fl is weak 
[see Figs. [5] (a), [5] (c), and [5] (e)]. With the interaction 
strength increasing [see Figs. E] (b), [5] (d), and [5] (f)], 
i.e., the dynamics of the Eq. © involving only the first 
and second-order terms in Eq. © are in good agreement 
with those obtained by the exact response (IC8I) . but the 
results obtained by the first-order approximation (linear 
regime) have serious deviations from the exact one (IC8I) . 
This difference comes from the nonlinear terms, which 
are ignored in linear response theory. 

Examining the Markovian regime, e.g., A = 15T, we 
find that the results given by the linear response (contain¬ 
ing only the first term in Eq. ©) are in good agreement 



FIG. 5: (Color online) n e (t) as a function of time. The result 
of linear response is given in (a), (c), and (e), while (b), (d), 
and (f) are for two-order nonlinear response. The red solid 
line, blue dashed line, and black dash-dotted line denote the 
exact infinity-order response (lC8l) . the linear response given 
by Eq. © containing only the first term of Eq. m , and 
the nonlinear response obtained by Eq. © containing only 
the first two terms in Eq. © , respectively. The parameters 
chosen are A = 0.2r,fl = 0.0ir,a = 4F,A = 0.2r for (a), 
A = 0.4F for (c), A = -0.2F for (e). f 1 = 0.2r, A = 0.2F for 
(b), A = 0.4F for (d), A = -0.2r for (f). 


with those obtained by the exact response (IC8D when the 
the interaction H is weak [see Figs. E] (a), [G] (c), and E] 
(e)]. When the the interaction strength H becomes strong 
[see Figs. [G] (b), [6] (d), and[6](f)], the dynamics given by 
Eq. © involving only the first and second-order terms 
in Eq. © are in good agreement with those obtained 
by the exact expression (1C8D . However, the results ob¬ 
tained by the first-order approximation (linear regime) 
have serious deviations from those obtained by the exact 
expression (|C8[I . The same observation can be found in 
the non-Markovian regime. 
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FIG. 6: (Color online) Comparison of the linear response 
[(a), (c), and (e)] and two-order nonlinear response [(b), 
(d), and (f)] in the Markovian regime. The red solid line, 
blue dashed line, and black dash-dotted line denote the ex¬ 
act response (IC8II . the linear response obtained by Eq. © 
contains only the first term in Eq. < 0 , and the nonlinear 
response obtained by Eq. © contains only the first two 
terms in Eq. respectively. The parameters chosen are 
A = 15T, a = 4T, A = T, H = 0.01T for (a), A = 5T,H = 0.1T 
for (c), A = -1T,H = 0.01T for (e). A = lT,fI = 0.5T for 
(b), A = 5T,fi = 5T for (d), A = -1T,H = 0.5T for (f). 


Appendix D: The derivation of Eq. (1281) . 


We first calculate the zeroth-order Hall conductivity at 
finite-temperature by rewriting Eq. m as 


(' n{p)\vn I (m(p)\v v \n(p)) 

+ Cnm 


p,m^n 

_ {fn{p)\vy. I n(p)) ( n(p)\v v \m(p)) 

frio e nm 

In the limit |w/e nm | <siC 1, 


(Dl) 


1 


1 


hu> ± i 


"nm °nm 


±i - 


hkJ 


+ 0(uj 2 ). (D2) 


Substituting Eq. (1D‘2I) into Eq. HD 111 , we have crffi (w) = 
cti(w) + cr 2 (cu) with 

ai ^ = ~Y Y fn(p)[( n (p)\ v v \m{pj) (m{p)\v v \n(p)) 

p,m^n 

+ {n(p)\v v | m{p)) (m(p)\v^ \ n(p))]/e„ m , 

(D3) 


and 


_ f Jl 

M fn(p)[( n (p)\ v n \ m (P)) {™>{p)\v v \n(p)) 

p,m^n 

- (n(j>)\v v | m(p)) (m(p)\v^ \n(p))]/e 2 nm . 

(D4) 

Now we show that the first term a\[ui) vanishes. Set 
[i — x and v = y), we have v x = x], then 


(n(p)\v x | m(p)) 


= ^(n(p)\[H,x\ \m(p)} 

'l 

= Te n m(n(p)\x\m(p )), 


(D5) 


and thus 

[(n(p)\v x | m(p)) ( m(p)\v y \n(p)) + ( n(p)\v v \m{p}) 

'i 

■ (m(p)\v x \n(p}}] = -e nm [(n(p)\x \m{p )) 

• < m{p)\v y | n{p)) - ( n(p)\v v \m(p)} ( m(p)\x \n(p))]. 

( D6) 

The factors e nm cancel each other, noticing 
X) \m(p)) (m(p)\ = I, we have 

m 

M = Y fnip){n{p)\[v y ,x] | n(p)) = 0, (D7) 

p,m^n 


since the commutator [x,v y \ vanishes, aj^J = 02 . As for 
the second term (1D4I) , simple algebra yields 



p,m^n 


fmnlm[( m (p)\v» I n{p)) (n(p)\y v I m(p))] 

^nm 


(D8) 

where f mn = [f m (p ) - fn(p)}- Therefore the finite- 
temperature Hall conductance for open system is given 
by cr^ = crffi + HT<r$. Following the same procedure, 
we have 


Y 9(0)fmnRe[{rn(p) |u M |n(p)) 
p,m^n 

■{n(p)\v v \m(p))]/el m + -^-y ^ h(0)f mn 

p,m^n 

■ (m(p )|tV | n(p}) ( m(p)\v v \n(p)} / ^ nm ' 

(D9) 

In the Heisenberg picture, the velocity operator (p = 
x , y , z) is defined as 

v r = = + (DIO) 


where the summations with respect to a automatically 
assumed in the Einstein summation convention. With 


Eq. (ID 101) 


(m(p) |u M \n{p)) 


dd a 

dpn 


(m(p)\a a | n{p )), 


(Dll) 
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for to ^ n, it is easy to prove that in the two-band model, Collecting all together, we can obtain the finite- 

temperature Hall conductance Eq. (051) for open system. 

Im[(m(p)|<7 Q | n(p)) (n{p)\ap \m(p))] =me a/ 3 7 -y, 

, / (D12) 

Re[(m(p)\a a \n{p)) (n(p)\a 0 \m(p))] = - 
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